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SIGNIFICANCE  AND  EXPLANATION 


Many  optimization  problems  in  engineering,  operations  research,  economics, 
etc.,  can  be  formulated  as  mathematical  programs , involving  the  extremization  of 
an  objective  function  under  certain  restrictions.  In  deterministic  mathematical 
programs,  the  data  are  assumed  to  be  given  exactly.  However,  in  many  situations, 
it  is  essential  to  take  into  account  probabilistic  (stochastic)  aspects,  e.g.,  the 
information  about  the  data  may  be  given  in  the  form  of  probability  distributions. 

One  way  of  solving  such  stochastic  mathematical  programs  is  to  reduce  them 

to  deterministic  mathematical  programs  whose  solution  provides  the  solution  of 

the  original  problem  in  a probabilistic  sense.  The  standard  deterministic 

T 

linear  program  is  to  maximize  c x subject  to  x > 0,  Ax  < b.  This  paper  con- 
siders the  case  where  c is  replaced  by  a vector  £,  the  components  of  which 
have  given  mean  values  and  a given  variance-covariance  matrix.  It  is  shown  that 
a certain  deterministic  quadratic  programming  problem  leads  to  a unified  treatment 
of  several  deterministic  models  proposed  for  this  problem  in  the  literature.  This 
approach  also  leads  to  an  efficient  computational  method  similar  to  the  well-known 
simplex  method  for  all  of  these  deterministic  models. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  summary 
lies  with  MRC , and  not  with  the  authors  of  this  report. 


A UNIFIED  PARAMETRIC  QUADRATIC  PROGRAMMING  SOLUTION 
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TO  SOME  STOCHASTIC  LINEAR  PROGRAMMING  MODELS 

* 

Masakazu  Kojima  and  Ikuyo  Kaneko 

1 . Introduction. 

The  problem  we  shall  consider  in  this  paper  is: 

T 

PO  maximize  £ x subject  to  x e X , 

where  £ is  an  n-dimensional  multi-variate  normal  random  vector  with  mean  c and  variance- 
covariance  matrix  V,  and 

X - { x € Rn  : Ax  < b,  x > 0} 

with  an  m by  n constant  matrix  A and  an  m-dimensional  constant  vector  b.  All  through 
the  paper  we  assume: 

(1)  V is  positive  definite:  and 

(2)  the  following  linear  program  has  a finite  optimum1: 

T 

Pi  maximize  c x subject  to  x e X 

One  approach  to  PO  is  to  consider  a deterministic  programming  model  which  produces  an 
"optimal"  solution  to  PO  in  a certain  sense:  various  approaches  in  this  direction  have  been 
proposed  (see  below) . These  models  differ  from  one  another  in  the  concept  of  optimality  in  the 
stochastic  program  PO.  In  this  paper,  we  shall  show  that  the  mean-variance  (parametric  qua- 
dratic programming)  approach  (of  Markowitz  [7])  yields  solutions  to  several  other  determinis- 
tic models  of  PO,  in  a unified  way. 

The  significance  of  the  result  is  two-fold.  It  clarifies  the  relationships  among  the 
solutions  of  these  models;  it  also  provides  an  efficient,  unified  computational  method  to 

^e  could  relax  this  assumption  somewhat  by  allowing  the  case  where  X is  empty:  the  assumed 
feasibility  of  PI  (or  PO)  is  merely  for  convenience. 
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solve  PO  based  on  these  various  models. 


[ 


The  deterministic  models  for  PO  we  shall  deal  with  in  this  paper  are  summarized  as 
follows: 

2 

P2  Find  efficient  solutions  of 

maximize  f(x)  subject  to  x e X , 
where  f(x)  = (c  x,  - x Vx)  . 

T 1 T 

P3(a)  maximize  ac  x - — x Vx  subject  to  x £ X , 

where  a > 0. 

T 

P4 (8)  maximize  Pr{£  x > B)  subject  to  x £ X , 

where  8 £ (-'»,“>)  ; and 

T 

P5(y)  maximize  0 subject  to  Pr{£  x > 0}  > y and  x e X , 

where  ye  (0,1) . 

The  problem  P2  is  the  so-called  mean-variance  model  due  to  Markowitz  [7] : it  can  be 
shown,  under  the  assumptions  stated  above,  that  all  efficient  solutions  of  P2  are  generated 
by  P3(a)  for  all  a > 0 (see  Theorem  2 below).  Solving  P3(a),  parametrically  for  all 
a > 0,  can  be  done  in  an  efficient  manner,  by  using  one  of  several  existing  methods.  For 
our  purpose,  the  parametric  version  of  Graves*  principal  pivoting  algorithm  (Graves  [3);  see 
Cottle  [1]  for  the  parametric  version)  will  be  most  convenient.  It  is  also  noted  that  the 
model  due  to  Freund  (2)  based  on  a risk  averse,  exponential  utility  function  amounts  to  solv- 
ing P3(a)  for  a > 0. 

The  model  P4(B)  is  often  referred  to  as  Charnes-Cooper * s P-model.  It  is  not  difficult 
to  show  that  under  the  condition  (1)  above  and  that  X does  not  contain  the  origin,  P4(B) 
is  equivalent  to  the  fractional  program: 

T 

c x-ft 

P4* (B)  maximize  — — subject  to  x e X 

/ T 

/x  Vx 


2A  ooint  x e x is  efficient  if  and  only  if  there  is  no  x e X such  that  f(x)  > f(x)  and 
f (x)  * f (x)  . 

-2- 
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It  is  worthwhile,  and  interesting,  to  observe  that  a certain  portfolio  selection  problem 
based  on  Tobin's  risk-free  assets  assumption  (see  [9])  is  formulated  (Lintner  [6])  as  the  same 


fractional  program  stated  above. 

The  problem  P4(B)  for  fixed  B may  be  solved  by  using  some  nonlinear  programming 
algorithm  applied  to  P4 ' ( 8) , or  alternatively  to  a certain  (strictly  concave)  nonlinear  pro- 
gram (P4"(B,X),  as  stated  in  Section  2)  equivalent  to  P4'(B).  The  problem  P4'(B)  arising 
in  the  risk-free  portfolio  selection  problem  can  be  transformed  into  a linear  complementarity 
problem  (see  Ziemba  et  al.  [10]  and  also  Pang  [8])  with  a oositive  definite  matrix;  and  thus 
it  can  be  solved  efficiently  by  some  complementary  pivoting  algorithm.  The  validity  of  this 

transformation,  however,  depends  critically  on  the  assumption  that  the  problem  contains  the 
T 

constraint  e x = 1,  where  e is  the  summation  vector  (which  is  a reasonable  assumption  in 
the  portfolio  setting) . 

Finally,  the  model  P5(y)  was  proposed  by  Kataoka  [5]^,  who  outlined  two  algorithms 
which  are  potentially  capable  of  solving  P5(y)  for  fixed  y.  These  algorithms  are  iterative 
schema  solving  a (concave  maximization)  quadratic  program  in  each  iteration.  In  [5]  it  was 
assumed  (in  addition  to  (1)  and  (2)  above)  that  0 ^ X and  y e (-j,  1)  . 

In  the  next  section  we  shall  derive  some  properties  of  the  solutions  of  P3(a)  and  also 
relate  P3(a)  to  P2.  We  shall  then  demonstrate  that  solving  P3(a),  a > 0,  will  simulta- 
neously solve  P4(8)  (treated  in  Section  3)  and  P5(y)  (Section  4),  for  all  relevant  values 
of  8 and  y respectively.  It  will  also  be  shown  that  solving  P4(8)  (or  P5(y),  re- 
spectively) for  fixed  value  of  6 (or  y)  can  be  done  efficiently  without  having  to  solve 
P3(a)  parametrically  for  all  a > 0.  Certain  relationships  between  the  values  of  the  para- 
meters (B  and  y)  and  the  corresponding  optimal  objective  values  of  the  problems  P4(6) 
and  P5(y)  will  be  expounded  in  the  fifth  and  final  section. 

2.  Some  Properties  of  the  Solutions  of  the  Parametric  Quadratic  Program  P3(g)  and  Their 
Relationship  to  Efficient  Solutions  of  P2. 

We  first  note  that  it  follows  from  (1)  and  (2)  that  for  every  a > 0,  P3(a)  has  a 

3 T 

Kataoka's  model  has  P(5  x > 0}  = y rather  than  > as  in  P5(y);  but  these  two  are  equi- 
valent under  the  condition  0 | X assumed  in  [5] . 

-3- 
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unique  solution.  In  fact,  under  these  conditions,  X is  nonempty  and  the  objective  value  of 
P3(a)  is  bounded  from  above  by  the  optimal  objective  value  of  PI;  uniqueness  is  an  obvious 
consequence  of  the  assumed  positive  definiteness  of  V. 

For  each  a > 0,  we  shall  let  <^>( oc)  denote  the  unique  solution  of  P3(a).  It  is  noted 
that  0(a)  for  all  a > 0 may  be  computed  by  the  parametric  principal  pivoting  algorithm;  in 
particular,  0( a)  is  continuous  and  piecewise- linear  in  a > 0. 

In  the  remainder  of  the  paper  we  shall  use  the  following  notation: 

X1  = {x  : x solves  Pi)  , 

2 

X = {x  : x is  an  efficient  solution  of  P2)  ; 

X3  = { «#>(a)  : a > 0}  , 

* T 1 

x = argmin{x  Vx  : x e X } , 

* 

a = the  final  "critical  value"  generated  in  the  parametric 

* 

Graves'  algorithm:  if  no  pivot  is  necessary,  set  a =0. 

T IT 

r(a)  * c v’(a) «£(a)  V<^(a)  (used  in  Section  3)  , 

a i 

-1  T 2 

s(a)  = — ( ifi(a)  V*s>(ct))  (used  in  Section  4) 
a 

* i 2 * 

Note  that  x exists  and  is  unique,  and  that  X n X = {x  }. 

We  shall  prove: 

Theorem  1 . 

* * 3 * 

For  a > a , 0(a)  = x and  (hence)  X = { 0(a)  : 0 < a < a }. 

Proof:  We  first  show  that  X3  is  bounded.  For  every  a > 0,  we  have: 

T 1 T T * 1 * T * 

ac  ^(a)  - — *(a)  Vv’(a)  > ac  x - — (x  ) Vx  , 

and 

T T * 

ac  ip(a)  < ac  x 


Since  V is  positive  definite,  the  boundedness  follows. 

* 

Now,  by  definition  a is  the  last  critical  value  generated  in  the  parametric  version 
of  Graves'  principal  pivoting  algorithm  (see  Cottle  [11)  applied  to  the  parametric  linear 
complementarity  problem  representing  the  Kuhn-Tucker  conditions  of  the  parametric  quadratic 


program  P3(a),  a > 0.  That  is,  the  same  complementary  basis  of  (I,  -M)  is  feasible  for 


every  a > a , where 


f T 

M = V A 
^ -AO 


It  then  follows  that. 


= ^(a  ) , a > a , 


i.e.,  ■*>( a ) solves  P3(a)  for  a > a . 

* 

For  a > a , a > 0,  we  then  have  that 

m m * l * m * 

(5)  c Ifi(a  ) - — ^(a  ) V v( a ) > c x - — (x  ) Vx 

2a  = 2a 

T * t * * * 

Letting  a -*■  “ in  (5)  we  obtain  c ^(x  ) > c x . But  since  x solves  PI  and  <?(a  ) e X, 

T * T * * 1 

the  inequality  holds  in  the  other  direction  also,  implying  c f(a  ) = c x , or  <fi(a  ) e X . 


With  this  equality,  (5)  yields 


* 7*  * * *T • * 

^>(a  ) V^( a ) < (x  ) Vx 


Thus  by  the  definition  of  x , we  have  a ) = x . 


Using  Theorem  1 we  shall  give  a proof  of  the  well-known  fact  that  all  efficient  solutions 


of  P2  are  generated  by  solving  P3(a)  for  a > 0. 


Theorem  2 . 


X = { *>(a)  : 0 < a < a } . 

2 3 2 3 

Proof : In  view  of  Theorem  1,  we  only  need  to  show  X = X . The  relation  X :>  X is  ob- 

2 

vious.  To  show  the  other  implication  let  x € X . Define  the  closed,  convex  set  K by 


4 


It  follows  from  x e X^  that  K has  the  empty  intersection  with  the  interior  of  the  non- 

2 T 2 

negative  orthant  of  R . Then  by  the  separation  theorem,  there  exists  p = (p^,  p^)  e R , 

T 

p > o,  p * 0 for  which  p y < 0 for  every  y e K.  Since  for  each  x e X the  point 


“ T > 

C X - C X 

1 T »T  * 

— [-x  Vx  + x Vx] 

belongs  to  K,  it  follows  that  for  every  x e X, 

T-  1 -T  - TIT 

PjC  x - — P2x  VX  > p^c  x - PjX  Vx 

- * 1 * 2 

If  p2  > 0,  then  x = ^(p^/Pj)  ; otherwise  (p^  > 0,  necessarily)  x e X . Since  x e X 

„ * - 3 

by  hypothesis  x = x and  hence  x e X . 

IS 

Finally,  we  shall  present  a couple  of  properties  of  V’(a)  which  will  be  useful  in  the 
subsequent  sections. 

Property  3 ■ 

T 

(i)  <^(ct)  V^(a)  is  monotone  nondecreasing  in  a > 0. 

(ii)  Assume  0 e X.  Then  either  ^(a)  = 0 for  each  a > 0,  or  ^>(a)  * 0 for  each 
a > 0. 


Proof:  The  proof  of  (i)  is  contained  in  Kataoka  [5) . To  show  (ii)  assume  0 e X 

(so  ^(0)  = 0)  and  suppose  ^(a)  * 0 for  some  a > 0.  From  (i)  and  Theorem  1 it  follows 

* 

that  x ^ 0 , which  implies  that 

T 

max{c  x : x e X}  > 0 ; 

* * * 

thus  c x >0.  Since  both  0 and  x belong  to  X,  so  does  Ax  for  every  A e [0,11 . 

For  A r (0,1),  we  have  that  the  objective  value  of  P3(a)  for  each  a > 0 is 

-6- 
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T * A * T * 

Aac  x — (x  ) Vx 

e T * 

Since  c x > 0,  this  value  is  strictly  positive  for  sufficiently  small  A > 0.  Thus,  x = 0 
can  not  be  the  solution  of  P3(cx)  for  a > 0. 

IS 

3 . Relationships  between  the  Parametric  Quadratic  Program  P3(g)  and  P-Model  P4(8). 

To  simplify  the  argument,  we  shall  first  assume  that  X does  not  contain  the  origin  and 
derive  the  relationships  between  P3(a)  and  P4(3).  We  shall  then  show  how  to  modify  the 
argument  to  handle  the  case  0 e X. 

As  was  mentioned  in  Section  1,  the  problem  P4(8)  for  fixed  8 is  equivalent  to  P4'(B) 

(under  0 ^ X) . By  using  a standard  result  in  fractional  programming  (see  e.g.  Jagannathan 
[4]),  P4'(B)  is  further  transformed  into: 

i 

T T 2 

P4"(B,A)  maximize  c x - 8 - A(x  Vx)  subject  to  x e X 

More  specifically,  for  fixed  8 e.  R x e X solves  P4'(8)  if  and  only  if  there  exists  A e R 

1. 

T T 2 

such  that  x solves  P4"(S,X)  and  c x - 6 - X(x  Vx)  = 0. 

We  prove: 

r 1 

Theorem  4 . 

Assume  0 ^ X. 

T 

(i)  For  every  a > 0,  'P(ct)  solves  P4(B)  with  8 = r(a)  , where  r(a)  = c via)  - 
«0(a)TV^(a)  /a . 

* * T * 

(ii)  For  every  8 e (-°°»  8 ) , where  8 = c x , there  exists  a > 0 such  that  8 = r(a) 

and  >8(  a)  solves  P4(8). 

* 

(iii)  For  every  x solving  P4(8)  with  8 e (-«•»  8 ),  it  holds  that  x = v»(a)  , where 

a = (xTVx)/(cTx  - 8). 

Proof:  (i)  Let  a > 0 be  given.  Since  •f(a)  solves  the  quadratic  program  P3(a),  the 

following  Kuhn-Tucker  conditions  hold: 

T T 

(6.1)  z = -ac  + V<8( a)  + A u > 0,  ^(a)  >0,  z *>( a)  = 0 , 

-7- 
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(6.2) 


Letting  X = 


y = b -A>0(a)  >0,  u > 0,  y u = 0 


T 2 

-s(a)  = (‘0(a)  VP( a))  / a > 0,  z'  = z/a  and  u'  = u/a,  we  have: 


z*  = -c  + Av,( a)/(  ^(a)Tv^(a) ) 2 + aTu'  > 0 , 


Y>(a)  > 0,  ( z * ) «2>(a)  = 0 , 


y = b - Ay>(oO  >0,  u'  > 0,  yTu'  = 0 , 


which  implies  that  ifi(a)  satisfies  the  Kuhn-Tucker  conditions  for  P4"(8»A).  Since  X > 0, 

the  objective  function  of  P4"(B,A)  is  strictly  concave  and  hence  ifi(a)  solves  P4"(B,A). 
T T 

From  B = c v>(a)  - (<fi(a)  Vv>(a))/a  it  follows  that 

1 

c%(a)  - 6 - X(^(a)TWia))2  = 0 


and  so  *P(a)  solves  P4'  (3)  , and  hence  P4(6).  This  proves  the  first  part. 

(ii)  The  second  part  is  proved  by  examining  some  properties  of  r(a)  . By  Theorem  1, 

* * 

Y>(a)  = x for  a > a ; and  so 

* 

lim  r(a)  = 6 


Since  0 ^ X,  ^(0)  * 0;  and  thus 


lim  r(a)  = -■» 
a-+0+ 


Since  r(a)  is  clearly  continuous  we  have  that 

(7)  {r(a)  : a > 0}  = (-<»,  8*) 

* 

Thus,  for  every  8 £ (-<*,  8 ) there  exists  a > 0 such  that  6 = r(a)  : by  part  (i)  ^(ct) 
solves  P4(8) • This  completes  the  proof  of  part  (ii) . 

(iii)  The  proof  of  part  (iii)  can  be  done  by  tracing  the  steps  of  the  proof  for  part  (i) 
in  the  opposite  way.  That  is,  if  x solves  P4(B),  then  it  solves  P4"(B,A)  with 


Furthermore,  it  can  be  shown  (see  Section  5)  that  r(a)  is  monotone  nondecreasing  in  a > 0. 


I 


A = (cTx  - 8)/(xTVx)2  > 0.  Then  by  setting  a - (xTVx)  2/X  = (xTVx)/(cTx  - B)  it  follows 
that  x sc.-.es  P3(a)  since  P3(a)  has  a unique  solut-on. 

IS 

The  above  theorem  implies  that  the  set  of  solutio  is  of  P3(a),  a > 0,  is  exactly  the 

* 

same  as  that  of  P4(B)  , 6 £ B ).  CNote  that  from  definition  it  follows  that  when 

* * Ik 

8 = 8 , x = <^(ct  ) solves  P4(8).)  In  view  of  Theorem  2,  it  may  be  argued  that  solving 

★ 

P4(8)  for  8 > 8 is  unnecessary,  at  least  from  the  standpoint  of  efficient  solutions. 
The  following  result  is  proved  directly  from  the  above  theorem. 


Assume  0 

4 x- 

Let  B e 

★ 

B ) and  a > 0 be  given. 

If  6 = r (a) , 

then 

•fi(a) 

solves 

P4 CB)  - 

if 

J > r (a)  , 

then  there  exists  a e (a,”) 

such  that  B = 

r (a) 

and 

a) 

solves 

P4(£f)  . 

if 

B < rla), 

then  there  exists  a e (0,  a) 

for  which  8 = 

r (a) 

and 

•f(a) 

solves  P4 ( 8 ) - 

This  result  is  useful  when  one  wants  to  solve  P4(8)  for  fixed  value  8-8;  it  suggests 
that  one  does  not  have  to  solve  P3(a)  for  all  a > 0.  Instead,  first  make  a guess  on  a 
likely  value  of  a solution  of  the  problem  P4(8)*  say  x.  In  many  practical  cases,  a reason- 
ably good  estimate  of  an  optimal  solution  is  available.  Let  a be  defined  by 

— T — 

— x Vx 

ct  = — > 0 , 

T — — 

c x - 8 


and  solve  P3(a) . Depending  on  the  relative  magnitude  of  r(a)  and  &,  one  then  solves 
P3(a)  parametrically  for  a > a,  or  a < a.  If  the  initial  guess  x is  close  to  the 
correct  solution,  this  second  stage  of  the  computation  would  generate  a such  that  8 = r(ot) 
quickly . 

It  is  worth  noting  that  the  value  of  r(a)  for  each  a > 0 may  be  determined  without 

evaluating  the  quadratic  form  of  <p(a)  . In  fact,  it  is  well-known  that  if  (x,  u)  is  a 

T T T 

Kuhn-Tucker  point  of  the  quadratic  program  P3(a),  then  it  follows  that  x Vx  = c x - b u. 
Thus  if  we  let  u(a)  be  the  dual  variable  in  P3(a)  corresponding  to  ^(a) , we  have  that 
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X; 


r (a) 


1 T 

— b u(a) 


In  the  remainder  of  this  section  we  deal  with  the  case  where  the  origin  is  contained  in 
X;  i.e.  b > 0.  In  this  case,  some  of  the  arguments  employed  above  are  no  longer  valid:  for 
instance,  the  problem  P4'(8)  is  not  well-defined  at  x = 0,  and  r(a)  does  not  approach 
-00  as  a approaches  to  zero.  However,  it  will  be  shown  that  essentially  the  same  results 
as  those  in  Theorem  4 and  its  corollary  hold  in  this  case  also. 


Theorem  5. 


Assume  0 € X. 


(i) 

For 

every 

ct 

> 0 such  that  r(a)  ^0, 

a) 

solves 

P4(8)  with 

e = r (a)  . 

(ii) 

For 

every 

8 

€ ( -°° , 0]  , x = 0 solves 

P4<0) . 

For  ' 

every  8 € (0, 

★ 

8 ) / there  exists 

o 

V 

o 

such  that 

8 

= r(a)  and  <p(a)  solves 

P4(B) 

(iii) 

For 

every 

X 

solving  P4(8)  for  8 e 

(0,  6* 

) » it 

holds  that  x 

= via)  , where 

a = (xTVx)/(cTx  - 3)  . 

T 

Proof : (i)  Since  r(a)  & 0,  we  have  ^(a)  * 0.  We  note  that  since  r(a)  = b u(a)/a  and 

* T 

b > 0,  r(a)  ^ 0 implies  r(a)  > 0.  Let  e = e »p(a)  > 0,  where  e is  the  summation  vector. 

Then  ^(a.)  solves  P3(a)  with  the  additional  constraint 

(8)  eTx  > e , 

* 

for  every  e e [0,  e ].  From  Theorem  4 it  follows  that  v?(c t)  solves  P4(r(a))  with  the 

★ 

additional  constraint  (8) , for  each  e r (0,  e ] . This  implies  that  no  x e X,  x 0 is  a 
"better"  solution  than  <£(a)  in  P4(r(a)).  Since  r(a)  > 0,  the  objective  value  in 
P4(r(c0)  for  x = 0 is  obviously  zero  but  that  corresponding  to  via)  is  positive  since 
<fi(a)  * 0;  thus  is  a "better"  solution  than  x = 0 in  P4(r(a)).  Therefore,  ^(a) 

solves  P4(r (a) ) . 

(ii)  First  note  that  whenever  8 < 0,  x = 0 solves  P4(8)  by  definition  of  the  pro- 

* * T * 

blem.  In  particular,  if  x =0,  then  8 = c x =0  and  so  x = 0 is  the  trivial  solution. 

* * 

Henceforth  we  shall  assume  x 0.  To  prove  the  assertion  in  (ii)  for  8 e (0,  8 1 , we 

shall  examine  some  properties  of  r(a)  under  the  assumption  0 € X.  First,  the  continuity 

i I 
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of  r(a)  and  that  r(u)  approaches  0 as  a -♦  • are  true  whether  or  not  0 € X.  We  have, 
however , 

(9)  lim  r(a)  = 0 

a-»-0+ 

To  see  this,  recall  that  <£(«)  is  continuous,  piecewise- linear  in  a;  and  thus  for  small  a, 
^>(a)  takes  on  the  form  >p(a)  = ad  for  some  n-vector  d,  Then  we  have 

r (a)  = ctcTd  - ~ (a2dTVd)  = a(cTd  - dTVd)  ; 

* 

hence  (9)  holds.  Thus  for  every  8 e (0,  6 ),  there  exists  a > 0 for  which  6 = r(a);  by 
part  (i)  4(a)  solves  P4(8). 

(iii)  The  third  part  is  proved  in  a similar  way  to  that  of  the  proof  of  part  (iii)  in 

Theorem  4 for  the  0 ^ X case.  We  note,  however,  that  unless  X = {0},  x = 0 cannot  be  a 

* * 

solution  to  P4(B)  for  8 £ (0,  8 ) : if  X = {0},  then  8 =0  and  so  the  statement  is 
true  by  default.  For  x ^ 0,  solving  P4(6),  we  can  "cut  off"  the  origin  to  make  the  argu- 
ment based  on  the  Kuhn-Tucker  conditions  valid. 

B 


The  following  is  direct  from  Theorem  5. 
Corollary. 


Assume 

0 £ X. 

Let  6 

e (0, 

* 

8 ) and 

a 

> 0 be 

given.  If 

6 = r (a) 

, then  v>(a)  solves 

P4  ( 8) 

. If 

8 > r(a) 

, then 

there 

exists  a 

e 

(a , oo) 

such  that 

8 = r (a) 

and 

V>(a)  solves 

P4(8> 

. If 

6 < r (a) 

, then 

there 

exists  a 

e 

(0,  a) 

such  that 

8 = r (a) 

and 

14(a)  solves 

P4(8)  . 

4.  Relationships  between  the  Parametric  Quadratic  Program  P3(g)  and  Kataoka  Model  P5(y). 

Assuming  ( 1 ) — ( 2)  and  0 <j  X,  Kataoka  [5]  showed  that  P5(y)  for  fixed  y e (0,1)  is 
transformed  into  the  following  problem: 

1. 

T T 2 

P5'(S)  maximize  c x + 6(x  Vx)  subject  to  x e x 

Here  6 = I ^(y),  where  for  s c R 
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t 


P5(Y)  need  not  be  treated  explicitly  and  hence  we  shall  say  x (rather  than  (x,  '))  solves 

the  problem  P5(Y). 

We  have: 

Theorem  6. 

Assume  0 ^ X. 

T 

(i)  For  every  a > 0,  ■f(a)  solves  P5(y)  with  y = I(s(a)),  where  s(a)  = -(^(a)  v (a))/a. 

(ii)  For  every  y e (j,  1)  there  exists  a > o such  that  Y = I(s(“))  and  'f(a)  solves 
P5(Y) . 

(iii)  For  every  x solving  P5(Y),  Y £ 1),  it  holds  that  x = v’f01)  > where 

1 

a = (xTVx)2/I-1(Y) . 

Proof : (i)  By  definition,  V>(a)  satisfies  the  Kuhn-Tucker  conditions  (6)  for  P3(u). 

Letting  u ' = u/a , and  z ’ = z/a  we  have 

1 

z'  = - c - s(a)Vf(cO/(*’(ct)TV>/>(a))2  + ATu'  > 0 , 

*(a>  > 0,  (z')%(ct)  = 0 , 

T 

y = b - A ifi(a)  >0,  u'  > 0,  y u'  = 0 

That  is,  'fi(a)  satisfies  the  Kuhn-Tucker  conditions  for  P5 ' ( 6 ) for  6 = s(a).  Since 

s(a)  < 0 the  conditions  are  sufficient  for  optimality.  Thus  0(a)  solves  P 5 ' ( 6 ) and  hence 

it  solves  P5(y),  where  y = I(s(a)).  This  completes  the  proof  of  (i) . 

(ii)  To  prove  (ii) , we  shall  examine  some  properties  of  s(a)  . By  Theorem  1,  we  have 

* * 

V5( at)  = x for  a > a ; so, 

lim  s(a)  = 0 

a-*» 

5If  x contains  the  origin,  then  P5(y)  may  be  restated  as:  Find  x e x,  x * 0 such  that 
T — — 1 /jp  — ip  _ /T 

c x + I (y)/x  Vx  = max{c  x + I (y)Yx  Vx  : x e x}  > 0:  if  no  such  x exists,  then  x = 0 
solves  the  problem. 
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Since  4(0)  * 0, 


lim  s(a)  = -“  . 

a->-0+ 

Furthermore,  s(a)  is  clearly  continuous;  thus 

(s(a)  : a > 0}  = (-“>,  0)  , 


or, 

{ X (s (a) ) i j > o)  « (i  1)  . 

Thus  for  every  y t (y,  1),  there  exists  a > 0 for  which  y = I(s(a)):  by  part  (i),  4(a) 

solves  P5 (y ) . 

H 

Again,  it  may  be  argued  that  solving  P5(y)  for  y < y is  unnecessary  from  the  stand- 
point of  efficient  solutions.  We  note  that  for  y = -j,  any  solution  to  PI  solves  P5(y). 
We  note  that  letting  u(a)  denote  the  dual  variable  corresponding  to  4(a)  in  P3(a),  then 
we  have 

(10)  s (a)  = -j(cT4(a)  - bTu(a))2  . 


For  computational  convenience  (see  the  remarks  immediately  after  Corollary  to  Theorem  4 
in  the  preceding  section)  we  state: 

Corollary. 

Assume  0 ^ X.  Let  y f (j,  1)  and  a > 0 be  given.  If  y = I(s(a)),  then  4(a) 
solves  P5(y).  If  y < I(s(a)),  then  there  exists  a e (a,°°)  such  that  y = I ( s (at) ) and 
4(a)  solves  P5(y).  If  y > I(s(a)),  then  there  exists  a e (0,  a)  for  which  y = I ( s ( a) ) 
and  4( a)  solves  P5(y)  . 

Now  we  shall  consider  the  case  0 e X. 

Recall  that  14(a)  is  continuous  and  piecewise-linear  in  a > 0.  If  0 c X,  then 

4(0)  = 0 and  so  for  small  a > 0,  4(a)  = ad  holds  for  some  n-vector  d.  In  view  of  Proper- 

* 

ty  3,  d = 0 if  and  only  if  x =0.  Let 

1 

(11)  s a lim  s(a)  = -(dTVd)2 

a->-0+ 
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I(s  ) if  d + 0 


d = 0 . 


Theorem  7. 


Assume  0 e X. 

(i)  For  every  a > 0 such  that  s(a)  * 0,  v>(a)  solves  P5(y)  with  y = I (c ■;&)). 

* 1 * 

(ii)  For  every  y e [y  , 1) , x = 0 solves  P5(y)  . For  every  ye  (j,  y ) , there  exists 

a > 0 for  which  y = Its  (a))  and  ^(a)  solves  P5(y). 

(iii)  For  every  x solving  P5(y)  for  y e (|-,  y* ) , it  holds  that  x = y>(a)  , where 

I * 

T 2-1 

a = (x  Vx)  /I  (y)  . 

* T 

Proof:  (i)  Since  s(a)  # 0,  y>(a)  / 0.  Let  c = e *(a)  > 0.  Then  y>(ot)  solves  P3(y) 

with  the  additional  constraint 

(13)  eTx  > e 

* 

for  every  e e [0,  e ] . By  Theorem  6 *(a)  solves  P5(I(s(a)))  with  the  additional  con- 

* 

straint  (13)  for  every  e e (0,  e 1.  This  implies  that  no  x e X,  x * 0 is  a "better" 

solution  than  ^(a)  in  P5(I(s(a))).  The  objective  value  o.‘  x = 0 in  P5(I(s(a)))  is  zero; 

whereas,  that  of  ^(a)  is  r(a)  (c.f.  Section  5)  which  is  nonnegative  (c.f.  the  proof  of 

Theorem  5).  Thus,  y>(a)  is  at  least  as  “good"  as  x = 0.  Hence,  *(a)  solves  P5(I(s(a))). 

* * * 

(ii)  If  y = 1/2,  then  x = 0 and  6=0.  In  this  case,  for  any  x e X,  x * 0, 

T T 

the  mean  c x of  the  normal  random  variable  5 x is  nonpositive:  thus  the  objective  value 

in  P5(y)  corresponding  to  x e X,  x * 0 is  nonpositive.  Hence  x = 0 solves  P5(y)  tri- 

* 1 » 

vially.  Assume  y > — . Then  by  the  definition  of  y , d * 0.  It  follows  from  Property  3 

* * 

that  x * 0;  and  so  0 >0. 

Suppose  x = 0 does  not  solve  P5(y)  for  some  ye  (y  , 1) . Then  there  exists  x e X, 

x * 0 and  6 > 0 such  that  PrfS^  > 6}  > y.  Since  y > 6 < 6*.  Letting  6 = 6/2  we 

- * t — * — * 

have  0 < 0 < 6 and  Pr{£  x > 0}  > y > y . By  Corollary  to  Theorem  5 we  can  find 


4 


» * A A A A 

a e (0,  a ] such  that  3 = r(a)  and  y(a)  solves  P4(8).  It  follows  then  that 
T—  - . * 

PrlC  x > 8}  < I(s(a))  and  hence  y < I(s(ci)).  On  the  other  hand,  since  s(a)  is  nondecreas- 

* * 

ing  in  a (Property  9 in  Section  5)  and  by  definition  of  & we  have  that  y > I(s(a)) 

* 

for  every  a > 0;  thus  a contradiction.  Hence,  x = 0 solves  P5(y)  for  every  ye  [y  , 1)  . 
The  proof  of  (ii)  is  complete  if  we  note  (11)  - (12),  that  s(a)  is  continuous  and  that 
s(a)  •»  0 as  a approaches  infinity. 

(iii)  Part  (iii)  follows  from  a similar  argument  to  that  proving  (iii)  in  Theorem  6. 

KJ 

Corollary ■ 

— l * — — — — 

Assume  0 e X.  Let  y e (— , y ) and  a > 0 be  given.  If  y = I(s(a))  , then  <?( a) 

solves  P5(y).  If  y < I(s(a)),  then  there  exists  a e (a,®)  such  that  y = I(s(a))  and 

v>(a)  solves  P5(y).  If  y > I(s(a)),  then  there  exists  a e (0,  a)  for  which  y = I ( s ( ot ) ) 

and  ^(.a)  solves  P5(y) . 

5.  Relationships  between  P-Model  P4(B)  and  Kataoka  Model  P5(y). 

Comparing  the  definitions  of  P4(6)  and  P5(y)  one  would  note  that  the  roles  of  the 
parameter  and  optimal  value  are  reversed  in  the  two  problems,  suggesting  some  close  relation- 
ships between  them.  In  fact.  Theorems  4-7  imply  the  following. 

Property  8. 

T T T 

For  a > 0 let  r(a)  = c v>(cx)  - ^(a)  Ha) /a  and  s(a)  = -(^(a)  V^tall/a.  Then  for 
a > 0,  ^(a)  solves  P4(r(a))  and  the  corresponding  objective  value  is  I(s(a));  tf(a)  also 
solves  P5(I(s(a)))  and  the  corresponding  objective  value  is  r(a). 

If  we  note  that  the  function  I(s)  is  strictly  decreasing  in  s,  the  above  result 
implies  the  following  relationship  between  r(a)  and  s(g):  for  a^,  > 0, 

(14)  signjrfa^)  - r(a2>}  = sign{s(a1>  - s(a2>} 

In  closing  the  paper  we  shall  present  the  following  result  (which  was  used  to  prove  the 
second  part  of  Theorem  7) . 
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Property  9 ■ 

Both  r(a)  and  s(a)  are  monotone  nondecreasing  in  a > 0. 

* 

Proof:  If  x =0,  then  ^(a)  = 0 for  all  a > 0 and  so  r(a)  = 0 and  s(a)  =0  for 

* 

all  a > 0.  Henceforth,  we  assume  that  x # 0;  thus  v>(a)  * 0 for  a > 0.  We  shall 
establish  the  monotonicity  of  s(a);  then  by  (14)  the  monotonicity  of  r(a)  follows. 

Noting  that  s(a)  is  continuous  and  is  converging  to  zero  as  a -*■  “>,  the  monotonicity 
is  proved  if  one  can  show  the  following: 


(15) 


0 < “l  < a2 


s (ci^ ) = s(a2> 


implies 


s(a)  = s( Oj) 


a e [a^ , a2) 


To  show  this  assume  0 < < a2  and  s(oc^)  = s(a^) . Let  6 = s(a^).  Let  a c (a^,  a^) 


be  given;  define 


1 a,  - a. 


and 


X,  = - 
2 a 


2 “l 


Then  X^  > 0,  X2  > 0,  X^  + Xj  = 1 and  a = X^a^  + *2a2‘  Let  E be  the  minimum  e 

•p  * * 

and  e v>(a2>  : since  x * 0 by  hypothesis  and  cx^,  a 2 > 0,  e is  positive.  Using  the  equi- 


valence of  P3(a)  and  P5'(s(a))  where  0 ^ X,  we  have,  in  this  case,  that  ^(o^) 


and 


tf(a2>  solve  P5'(s(a))  with  the  additional  constraint 


T .. 

e x > e 


(16) 

* 

for  every  £ e [0,  e ) . 

Since  6 < 0,  P5'(6)  is  a concave  maximization;  thus  x = X^^(a^)  + X2^(a2>  also 


solves  P5'(S).  Then, 


T—  — t - 2 

c x + 5(x  Vx) 


l Xi(c%(ai)  + 6(  V>(cti)TV^(ai) ) 2 ) 

i=l 


= c x + l X^f^fcu)  v#( a/)) 
i=l 
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It  follows  then 


(i’vx)2  - l A (*(a  )TV*(a  ))2  . 

* * * 


By  definition  we  have 

1 

« - stdj)  - - (^(a.)TV^>(a1))2  , 


tl8>  (*<a.)TV*(a.))2  - - a. i,  i - 1,  2 . 

Substituting  the  quadratic  forms  in  (17)  by  (18)  we  obtain 

l 

(19)  a - s(o,)  - - - (x'Vx)2  . 

l a 

The  implication  (15)  will  be  proved,  then,  by  showing  that  x - <p( a)  . 

Since  x solves  P5 • ( a ) with  (16),  the  following  Kuhn-Tucker  conditions  hold  for  some 
u,  u , y and  y : 

1 

(20.1)  z = - c + avx/lii^Vx)2  + hu  - euQ  > 0 

(20.2)  x > 0,  zTx  = 0 , 

(20.3)  y = b - Ax  > 0,  u > 0,  yTu  = 0 , 

» T 

(20.4)  yQ  • - e + e x > 0,  u0  > 0,  yQ"uo  = 0 • 

1 

-T  — 2 

Multiplying  (20.1)  by  a (<=  - (x  Vx)  /6  by  (19)),  and  letting  u'  = au  and  u^  = auQ,  we 
see  that  (x,  u',  u^)  satisfies  the  Kuhn-Tucker  conditions  for  P3(a)  with  (16).  Letting 
s -►  0 we  conclude  that  x solves  P3(a);  since  P3(a)  has  a unique  solution  x = V>(a)  , 
as  needed. 


L 
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